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ABSTRACT 

We propose convergent plasma lenses, possibly from current sheets, as a generic solu- 
tion to strong interstellar scattering. These lenses resolve the overpressure problem by 
geometric alignment as noted by Goldreich and Shridhar (2006). They further quanti- 
tatively explain properties of extreme scattering events, and pulsar parabolic arcs. This 
model makes quantitative predictions testable by VLBI on scattering events. It differs 
conceptually from previous models by occurring through rare, localized underdense 
sheets. Such sheets are thermally and kinematically stable, and could be consequences 
of reconnection. The apparent diffractive effects are a result of coherent interference 
of refractive images. We propose that these lenses can be used for precision distance 
determination to pulsars, enabling accurate gravity source localization. 

Key words: Interstellar Medium, reconnection, extreme scattering events 



1 INTRODUCTION 

Several scintillation phenomena in the interstellar medium 
have posed challen ges to any physical model to explain 
them (|Rickett| l2007Y These include: 

(1) extreme scattering events (ESE) (jPiedler et al.l llQSTi ). 
Compact radio sources are occassionally observed to go 
through a period of demagnification at low frequencies by 

roughly a factor of 2; 

• J^ (2) pulsar parabolic arcs (|Stinebring et al1l200ll ): 
(3) galactic center scattering. 
5^ In each of these cases, a simple application of Snell's 

law with the assumption of spherical symmetry of the lens 
requires free electron densities up to ~ 10'* cm~^. Free elec- 
trons are at temperatures of at least ~ 10 K, and the in- 
ferred pressures are difficult to reconcile with pressure bal- 
ance in the interstellar medium. 

A solution to case (3) has been proposed by 
iGoldreich fc SridhaJ l|2006l ). who pointed out that scatter- 
ing for sheet-like structures is dominated by the ones most 
aligned with the line of sight. The alignment lowers the re- 
quired three dimensional electron density. 

In this paper, we compute the quantitative conse- 
quences of plasma lenses, and show that triaxial structures 
are consistent with all observational data without requiring 
any unusual physical conditions. Large axis ratios are generic 
consequences of reconnection. In ideal resistive MHD with 
ohmic conversion of magnetic fields, current sheets would 
be overdense in pressure equilibrium. Since the resistivity is 
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almost certainly not ohmic, the actual density is not known. 
The phenomelogy suggests underdense current sheets, which 
is a probe of the physics of reconnection. Geometric factors 
cause the scattering to be dominated by rare aligned events. 



2 PLASMA LENSES 



We follow the notation of lSchneider et al] l|l992l ). reproduc- 
ing their lensing geometry in Fig.[T] The diameter distances 
from the observer to the lens plane and to the source plane 
are Dd and Da respectively, and the distance of the source 
plane from the lens plane is Dds- Physical coordinates in the 
source and lens planes are ff and £, respectively, defined with 
respect to the optic axis connecting the observer with the 
centre of the lens. The deflection of a light ray at the lens 
plane is denoted by a. The angular position of a source at 
ff is l3, and 6 is the apparent angular position from which 
the deflected ray travels. They are related through the lens 
equation: 



/3 = e- 5^d(Dd0) = e- a{e) 

Us 



(1) 



where a is the deflection angle, simply the difference between 
the angular position of the source and image. 

The simplest geometric lens we consider is a triaxial 
G aussian electron dens ity distribution. This is an extension 
of'Roman i et all (119871 ) , where we consider the limit that the 
axis ratios are very large, corresponding to thin sheets. The 
Romani et al picture was further quantified bv lClegg et al.l 
l|l998l ). The main qualitative difference in our new analysis 
is to consider the convergent case, which solves several key 
problems. 
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Figure 1. Geometry of a lensing event (reproduced from Schnei- 
der et al. 1992). See text for details. 



In projection, sheets are highly elongated two dimen- 
sional Gaussian surface densities. The lensing physics is thus 
two-dimensional with ^ chosen to be aligned with the short 
axis, with the deflection potential tj) (the projection of the 
3-D Newtonian potential) given by: 



il){9) = (TgKo exp(-6'^/2crg) . 



(2) 



where the gradient of the potential, Ve "ip, gives the deflec- 
tion angle. For a convergent lens, kq < 0. 

The phase velocity and group velocity of radio waves in 
a plasma are 



Cph = C/ ^/(T 

and 



Cg = cy/ (1 -a;2/a;2) 
where 



(3) 



(4) 



cjp = ^neC^/eome (5) 

with electron number density n^, electron charge and mass 
e and rrie respectively, and eo being the permittivity of free 
space. Note that CphCg = c^. In the ISM, typical plasma 
frequencies are kilohertz, so the frequencies of relevance are 
much larger than the plasma frequency, w ^ Wp. 

The effective refractive index of the plasma and the 
potential are related through n = 1 — 2$/c^ = c/cph, 
and thus the intrinsic component of path (or phase) delay, 
due to passage through the plasma is frequency-dependent 
CTgrav = / {n{uj) — 1) dz. This refractive delay is the anal- 
ogy to the gravitational Shapiro delay. In the limit lu ^ uup, 
use the notation of gravitational lensing. 



(6) 



where the first term accounts for the geometrical delay due 
to the offset in source and image positions. 

Two possibilities exist: electron overdensities result in 
a faster phase velocity, corresponding to a c oncave (diver- 
gent) optical lens. This case was considered bv lRomani et al.) 
([1983), which results in 4 sets of caustics, which are not con- 
sistent with the observed properties of ESEs. At the center of 
one of these events, only one image exists, with flux oc 
while observed events have only order unity flux reduction. 
Physically, electron overdensities can also lead to cooling 
instabilities, and are difficult to pressure confine. 

Electron underdensities are convergent lenses. These are 
generic consequences of heating processes, for example mag- 
netic reconnection events. At the centre of these lenses, as we 
will demonstrate, three images can be seen: a faint central 
image and two brighter ones on each side of the lens. 

The convergence of the lens is related to the laplacian of 
the deflection potential, i.e. k = tp" /2 and the components 
of the shear for a sheet with the short axis on the plane of 
the sky are 71 = k, 72 = 0. Thus the magnification fi — 
1/(1 — 2k). The convergence is a dimensionless measure of 
the distance to the lens in units of its focal lens. When the 
convergence is small, the lensing is weak, and only one image 
is formed. 

The critical points in the lens plane occur when 1/fj, — 0. 
The previous treatment by Cordcs ct al. (1986) neglected 
the fact that the source plane position and the lens plane 
positions differ, leading to qualitatively different behaviour 
in the strong lensing regime. In analogy with gravitational 
lensing, we call the general lensing geometry "strong", to 
distinguish from the "weak" lensing limit which is easier to 
compute and has been primarily used in the scintillation 
literature. In weak lensing, the image plane and lens plane 
are equivalent, while strong lensing considers the general 
case. 

For the Gaussian lens, we can solve for the critical 
points and caustics in the limit of large convergence: 



± 1 + 



2ko 



I3c = ± 



1^0 



(7) 



Note that here e ~ 2.718 refers to the natural log base, not 
the electron charge as before. 

On axis, the three images for /3 — appear a.t 6 = 
0, ±2 log(— Ko)- Their magnifications are 

-1 1 



1 + Ko ' 21og(— Ko) — 1 ■ 



(8) 



Note that a convergent lens has three images: two bright 
and one faint. These could be observed with VLBI, or 
through pulsar scintillation. 

We show a light curve in Figure [21 for a source moving 
parallel with respect to the short axis of the sheet. As a 
function of source position, or time with respect to direct 
alignment between source and lens, the lower panel shows 
the total image magnification for a source which is half the 
size of the lens, the middle panel shows the magnification of 
the individual images for a point source, and the upper panel 
shows the corresponding positions of the images (indicated 
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Figure 2. Multiple image creation and magnification during an 
extreme scattering event. As a function of source position or time, 
top panel shows positions of image(s), middle and bottom panels 
show their individual and total magnifications respectively. The 
bottom panel is for a Gaussian source profile of half the angular 
size of the lens. Negative magnifications correspond to an image 
parity change, which becomes observable in terms of the pulse 
rotational reflex motion. 

with the same line styles). Note the production of 3 images 
and boost in magnification upon caustic crossing, and the 
net total demagnification during a large period of the event. 
This is much like looking at a magnifying glass at distances 
longer than the focal length: the central image is inverted, 
and demagnified. 

Along a path with many randomly oriented lenses, the 
convergence is dominated by the lenses where the interme- 
diate axis is aligned to the line of sight, and in projection 
appear of comparable size to the short axis. Since the lenses 
are rare, a general line of sight may not contain any such 
highly aligned lenses. Neverthe less, the lensing is dominat ed 
by the most aligned scatterer l|Goldreich fc Sridharlliooil ). 



3 PULSARS 

The multiple images of the lens can be thought of as an 
interferometer on the source. Pulsars generally remain un- 
resolved under such lensing events, and the images interfere 
with each other. 

The time delay through a lens is t cc ip + o? jl. For a 
convergent plasma lens, i/; < 0, resulting in a negative time 
delay on axis, i.e. an advanced pulsar arrival. 

Fig-E] shows the offset in arrival-time delay of pulses 
from a pulsar during a lensing event. Note the sharp rise in 
delay near caustic crossing, and the apparent negative delay 
around the time when the pulsar and lens are in alignment. 
While the phase velocity in a plasma is larger than the speed 
of light, the group velocity is always smaller. In a convergent 
plasma lens, the underdensity results in a lower phase veloc- 
ity, and thus higher group velocity, resulting in the apparent 
advancement in the pulse arrival. 
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Figure 3. Arrival-time delay of pulses from a pulsar averaged 
over all apparent images. 




doppler velocity 

Figure 4. Pulsar image trajectory in the secondary spectrum. 
The triangle marks the position at caustic crossing, with one tri- 
angle superimposed under the central square. The squares mark 
the positions a short while later, 6% into the ESE. The central 
image has barely moved. The octagons mark the image positions 
at the center of the event. In a snapshot, only can only measure 
relative delays between images. With pulsar timing, the actual 
offset is also observable, leading to a negative group delay for the 
central demagnified image. 



3.1 Incidence of lensing 

The incidence of ESE is observed to be roughly 1% at 2.7 
GHz. At 400 MHz, this gives about a 2% probability of en- 
countering a lens with convergence >45. Each lens is visible 
at 45 radii, implying that half the time a pulsar dynamic 
spectrum can show an inverted parabola. 
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3.2 Distance measurements 

ESE lensing events in pulsars open up the possibility of pre- 
cise geometric distance determinations, which in turn can 
increase the sensitivity of gravitational wave detection, and 
result in high angular resolution (,Boyle_& Pen 2010). It re- 
quires two lenses at distinct distances. During an ESE, two 
bright images will have known mutual time delays, readily 
measurable to a part in 1000 or better. There is a reason- 
able chance of an offset scattering lens at a different dis- 
tance: these inverted parabolae are observed in a substan- 
tial fraction of all pulsars. Using interstellar holography, the 
resulting fringe pattern can be decomposed into four dis- 
tinct images: two from the ESE lens, and two more from the 
second scatter. VLBI allows a measurement of the angular 
separation of the two ESE images, and the angular distance 
to the second scatterer. These separations will be compa- 
rable to VLBI resolution, yielding fractional distance errors 
of ~ 1/SN, which can readily exceed 10'^ for bright pul- 
sars relevant to gravitational wave timing experiments. The 
effective distance to the ESE is measured. The two images 
will have different effective distances to the second scatterer. 
The three distances combined result in a measurement of the 
three unknowns: distance to the two lenses, and the distance 
to the pulsar. 

An additional option exists for distance to pulsars in 
binary systems. The pulsar performs a periodic orbit which 
is an ellipse on the sky. The known orbital parameters, and 
the time resolved lens measurements, allow a reconstruction 
of lens and sources distances. 



4 EXTRAGALACTIC RADIO SOURCES 

The original ESE were found in extragalactic radio sources. 
The classic example is 0954-1-658. 

At 8.1 GHz, a complex light curve containing multiple 
peaks was observed. At 2.7 GHz, the light curve much re- 
sembles the bottom panel of Figure (2] In our picture, we 
consider that kq ~ 1 at 8.1 GHz. For frequencies close to 
the critical, small substructures in the lens or the source 
can lead to rapid flux variations, which is observed. At 2.7 
GHz, one then expects hio ~ 10, which was the value chosen 
in the figure. We see that the model matches the data well. 
The factor of two demagnification is a generic consequence 
of convergent lenses. Divergent lenses, on the other hand, 
result in a ~ 1/ko demagnification (see Eq[Sl), which would 
be an order of magnitude scaling from the near critical mag- 
nfication near 8.1 GHz. 

The parallactic alignment of the lens with respect to the 
source is not known, making it challenging to di rectly inter- 

J iret t he existing VLBI structures of the source. ILazio et al] 
2000l ) noted an apparent increase in angular size during the 
ESE, which they interpreted as inconsistent with refracive 
lensing, which would conserve surface brightness and require 
a decrease in angular size. In the Gaussian model, the ESE 
is associated with multiple images with separations smaller 
than the VLBI beam. Thus the apparent size increase is in 
fact a generic consequence of this model. 



5 GALACTIC CENTER 

The radio source Sgr A* and neighboring masers 
are observed to be substantially scatter broadened 
(jLazio fc Corded 1 19981 ). 6 ~ l"(GHz/i/)^ Models place the 
screen within ~ 150 pc of the galactic center, implying a 
scattering angle of order an arcminute. This has led to chal- 
lenges for confining enough plasma to create this scattering. 
iGoldreich fc Sridhaij (|2006l ) proposed sheets as a solution. 
Statistically, the scattering is always dominated by the small 
fraction of sheets that are seen face-on. The convergent lens 
picture works in this regime as well. The recently discovered 
pulsar J1746-2850II (|Deneva et al.ll2009l ') might be behind a 
lensing screen, and might exhibit frequency dependent lens- 
ing. Frequent monitoring would be required. VLBI measure- 
ments of the scattering disk could also yie ld insights. Th e 
scattering disk is elongated by about 2:1 (|Lo et al.lll993l l. 
This suggests that the number of scatters is either small, or 
the magnetic fields have substantial systematic alignment. 
Magnetic field lines must lie in the plane of the plasma sheets 
to prevent diffusion and thus lens broadening. Images are al- 
ways scattered perpendicular to the field line direction. 



6 OBSERVATIONAL TESTS 

Two regimes can be searched to confirm this model of these 
convergent lenses. The most direct is the ob servation of 
an ESE with VLBI. It had been done once by iLazio et al.l 
(2000), but with insufficient resolution to discern the mul- 
tiple images. They note a slight broadening of the source 
at the lensing peak, which might be due to an unresolved 
mixing of two images. Lensing conserves surface brightness, 
so one expects a larger image to be brighter. The authors 
used this argument to rule out a refractive divergent lens. 
We note, however, that the expected lens image separation 
is only modestly, say a factor of 2 larger than the source. 
This results in two images, each of which is demagnified 
and smaller, but next to each other to make them appear 
more extended. Their synthesized beam is much larger than 
the source size or the separation, so these two effects are not 
separable. 

Unfortunately, the observational monitoring pro- 
gramme of bright compact sources has been discontinued. 

Pulsars also u ndergo extreme scattering events 
jCognard et al.lll993l ). In such cases, the three images will 
coherently interfere with each other, resulting in a distinc- 
tive signature in the secondary spectrum, shown in Figure 
|4] This figure plots the raw delays which can be obtained 
from interstellar holographv (|Walker et al.]|2008l ). 



7 PHYSICAL DISCUSSION 

So far, the discussion has focussed on the quantitative prop- 
erties of the convergent lens. Some comments on its relation 
to the physical properties of the lenses are in order. 

Typical electron densities as measured by pulsar dis- 
persion measure are rie ~ 0.03. This gives a phase velocity 
about 10~^^ larger than c at frequencies of ~ GHz. For an 
underdense lens, the maximal deflection angle from Snell's 
law is therefore 10~*/a mas, where a is the angle between 
the lens surface and the line of sight. Grazing incidence with 
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a ~ 10~* is needed to explain the observed refractive angles 
of ~ mas. 

We speculate that reconnection current sheets could 
satisfy the properties required by the phenomenology. Un- 
fortunately, the quantitative process of reconnection is very 
poorly understood. The thickness of a current sheet is prob- 
ably related to resistivity. The primary quantitative calcu- 
lable ohmic restivity leads to very long reconnection time 
scales, leading to problems with the generation of magnetic 
fields. Alternative reconnection processes, or fast reconnec- 
tion, has been proposed. The reconnection rate and sheet 
geometry differs by many orders o f magnitude between pro- 
posed scenarios (|Pang et al.ll2010l '). Without a known recon- 
nection rate, the incidence rate of current sheets is not pre- 
dictable, nor is the lens thickness. Similarly, the aspect ra- 
tio of the sheets is not known. One qualitatively expects 
the longest axis to be related to the curvature of field lines, 
which determines the amount of area that approaches suf- 
ficiently closely to undergo reconnection. The shortest axis 
is given by the unknown resistive scale, which could be as 
short as the electron gyromagnetic radius. To explain the 
phenomenology of ESE, we require the ratio of the interme- 
diate axis to the short axis to be > 10'*. Unfortunately, the 
statistical properties of ESEs is not predictable within our 
understanding of plasma physics. Instead, we can use the 
properties of the ESE to learn about reconnection. 



8 CONCLUSIONS 

We have introduced triaxial convergent plasma lenses to 
explain three separate, but related, refractive plasma lens- 
ing effects. Reconnection events generically lead to highly 
fiattened current sheets. We have found that generic phys- 
ical conditions can account for all the known phenomena, 
which are otherwise challenging to explain. These include 
the apparent distance localization of scatters, the modest 
frequency dependence of ESE, overpressure problem. It also 
suggests that the apparent diffractive phenomena in pulsar 
may be due to interference of refractive images, rather than 
stochastic fresnel scale lensing. 
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